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Quantum thermodynamics can be understood as a resource theory, whereby thermal states are free 
and the only allowed operations are unitary transformations commuting with the total Hamiltonian 
of the system. Previous literature on the subject has just focused on transformations between 
different state resources, overlooking the fact that quantum operations which do not commute with 
the total energy also constitute a potentially valuable resource. In this Letter, given a number of 
non-thermal quantum channels, we study the problem of how to integrate them in a thermal engine 
so as to distill a maximum amount of work. We find that, in the limit of asymptotically many uses 
of each channel, the distillable work is an additive function of the considered channels, computable 
for both finite dimensional quantum operations and bosonic channels. We apply our results to 
bound the amount of distillable work due to the natural non-thermal processes postulated in the 
Ghirardi-Rimini-Weber (GRW) collapse model. We find that, although GRW theory predicts the 
possibility to extract work from the vacuum at no cost, the power which a collapse engine could in 
principle generate is extremely low. 


The field of quantum thermodynamics has seen a surge 
in interest in the past years, with increasing attention to¬ 
wards testing the validity of the rules of classical thermo¬ 
dynamics in the quantum regime. A major topic within 
thermodynamics is that of extracting work out of a given 
system and the optimal way to perform this. One way 
this has been approached in the quantum case was by 
considering it from the perspective of a resource theory. 

The idea of a resource theory of thermodynamics is to 
assume one has unlimited access to thermal baths (i.e. 
Gibbs states of a fixed temperature T), the freedom to 
apply any energy conserving unitary on system plus the 
bath (any unitary V that commutes with the total Hamil¬ 
tonian), and the possibility of discarding part of the sys¬ 
tem or bath (i.e. apply partial traces). These rules are 
imported from classical thermodynamics, where one as¬ 
sumes access to infinite baths of constant temperature 
and any evolution where energy is conserved. 

As a matter of fact, resource theories have been very 
useful in different topics within quantum information the¬ 
ory EHl. The idea is similar to above: considering 
free access to certain operations and/or states, any state 
and/or operation that is not in the above set can in prin¬ 
ciple be used as a resource. 

This work complements previous research in quantum 
thermodynamics by accommodating the possibility of 
considering non-thermal maps, or channels, as a resource. 

Physical operations are represented by quantum chan¬ 
nels, i.e., completely positive trace preserving maps Q : 
B['H) —>• B{H') acting on a state space B['H). For sim¬ 
plicity, we will assume that the input and output spaces 
(and, as we will see later, Hamiltonians) of each channel 
are the same, although the results can be easily general¬ 
ized. 

Unitary evolution is a particular instance of a quantum 
channel, determined by the evolution operator. However, 
quantum channels allow to express more general evolu¬ 


tions. For instance, any unitary interaction of a system 
with an ancilla (or environment) generates a quantum 
channel, given by 

Q(p) =Tr^[Up®cTAUt], (1) 

where a a is the state of the ancilla, and V is some unitary 
operator In fact, it can be shown that any channel can 
be generated via the above procedure . 

If, in the above expression, the ancilla is in a Gibbs 
state of temperature T and the unitary V commutes with 
the total Hamiltonian Hp = Hs ® 1a + Is <8) Ha of 
the target-ancilla system, the resulting map is called a 
thermal channel. These will be the free operations in our 
theory, while any non-thermal map Q will be considered 
as a resource. 

In this scenario, we define quantum work W as the 
process of exciting a two-level system with Hamiltonian 
H = IU|1)(1| from its ground state |0) to the excited state 
11) Different authors have explored how much work 
one can extract from a non-thermal quantum state i- 
[l^ . When we regard the maximum average work as a 
figure of merit, a quantum generalization of the classical 
free energy naturally emerges SSEl: 

F{p) = U{p) - KbTS{p). (2) 

Here p is the state of the system from which we wish 
to extract work; U{p) = TrjpTJ], its average energy; and 
S{p) = — Trjplogp], its von Neumann entropy. Kb is 
Boltzmann’s constant. The maximum amount of work 
one can extract (on average) from the state p can then 
be shown to be F{p) — F{Tth), where Tth represents the 
thermal state at temperature T. 

In this Letter we want to address the following re¬ 
lated problem: suppose we want to build a thermal en¬ 
gine, where we are allowed to integrate a number of non- 
thermal gates each of which is assumed to act 
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on a system with Hamiltonian Hi. More specifically, our 
machine can make free use of any amount of thermal 
states and operations, and we can invoke one use of each 
of the channels in any order we want at any step. 

We are also allowed to use catalysts, i.e., we can use any 
number of non-thermal states, as long as we return them 
in the end. Under these conditions, what is the maximum 
amount of work that our device can extract? 

There are two ways to approach this problem: 

1. We can restrict to thermodynamical processes 
which distill work deterministically, i.e., always the 
same amount. The corresponding deterministic ex¬ 
tractable work can then be shown to behave very 
badly: not only is it not additive, but it can be 
super-activated. That is, there exist channels H 
such that no work can be distilled from a single 
use, but two uses of the channel can be combined 
to produce a non-zero amount of deterministic work 
(see the Supplemental Information). 

2. Alternatively, we can consider thermodynamical 
processes which generate a given amount of work 
with high probability. Here the figure of merit 
would be the maximum amount of work that can 
be distilled asymptotically (on average) when we 
have access to n uses of each channel. 

We will follow the second approach: in the next pages 
we will show that the asymptotically extractable work is 
upper bounded by Hi), where 


W{n,H) =maxAF{p,n), (3) 

p 

with AF{p, H) denoting the free energy difference be¬ 
tween the states H(p) and p, i.e. 

AF{p,n) =TTmip) - p)H] - KsT[Sinip)) - S{p)]. 

( 4 ) 

The quantity W{Vl,H) will be called the distillable work 
of channel H. From the inequality AF^xth,^) > 0, it 
follows that W{Vt,H) > 0 for any fl. 

The bound W{Vli,Hi) can be achieved asymptot¬ 
ically via a simple protocol where we prepare suitable 
initial states Ucat (the catalysts) maximizing eq. ([3|) for 
each channel, and then let each channel act over its corre¬ 
sponding maximizer. The result of this protocol will be a 
state with free energy F{acat) + W{V,i, Hi). Given 
access to n uses of each channel, we can thus prepare n 
copies of the latter state, whose free energy can be con¬ 
verted to work via thermal operations using the protocol 
depicted in Q. Following Q, part of this work (roughly 
nF{acat)) can then be used to regenerate the catalysts 
up to a small error [l^ . The average work extracted with 
this procedure (namely, the total work divided by n) is 
thus given by Y)f=iW{^i, Hi). 


Note, though, that, unless the catalysts are diago¬ 
nal in the energy basis, an extra amount of coherence, 
sub-linear in n, may be needed to rebuild them (see 
Appendix E of 0 )- More specifically, for each energy 
transition ^ Et in the Hamiltonian Hi, the proto¬ 
col proposed in Q requires a system with Hamiltonian 
= Y2^\Es-Et)k\k){k\ instate ^ ELoI^)’ with 
m sublinear in the number n of uses of each channel. Like 
the catalyst states, at the end of the protocol such ‘coher¬ 
ent states’ will be approximately rebuilt with vanishing 
error. 

In order to prove the above result, and some later ones, 
the next lemma will be invoked extensively: 

Lemma 1. Let be an N-partite quantum state, and 

let be a collection of N single-site quantum chan¬ 
nels. Defining Hi. jv = we have that 


N 

E 

2 = 1 


5(a,)-5(U,(a,)) > 5(aW)-5(Hi,..^(a^"0). (5) 




The proof is a straightforward application of the contrac- 
tivity of the relative entropy 14| . 

An almost immediate consequence of Lemma [T] is that 
W (H, H), as defined by eq. ([3]), has the remarkable prop¬ 
erty of being additive. That is, if the bipartite system 12 
is described by the Hamiltonian H 12 = l 2 + li®H 2 , 
and the channels Hi and H 2 act on the respective Hilbert 
spaces then, lU(Hi (g) D. 2 ,Hi 2 ) = W{D.i,Hi) -\- 

W{Ll2,H2). 

Indeed, let H 12 = Hi 0 H 2 act on the bipartite state 
P12. By choosing pi 2 = pi® P2 in m we trivially have 
that VF(Hi 0 H 2 ,Lfi 2 ) > lU(Hi, iLi)-I-1U(H2, ^ 2 ), since 
maximizing over states in 12 is more general than maxi¬ 
mizing over 1 and 2 independently. Let us then focus on 
the opposite inequality. By Lemma [T1 we have 

^ ^(p,) - ^(H,(p,)) > ^(P12) - 5(Hi2(pi2)). (6) 


Substituting into ([3]) gives 

E/ ^E{pi,Lli) > AF(pi2,Hi2) V P 12 , (7) 

i=l,2 


It follows that J2i=i,2^i^^’Hi) > 1F(Hi 2, i7i2)- 

We are now ready to prove that W{fl,H) quantifies 
the maximum (average) amount of work one can extract 
from channel H. 

Proposition 1. Let be a set of quantum chan¬ 

nels, defined over different quantum systems with Hamil¬ 
tonians {Hi}fLi. Suppose that we integrate n uses of all 
such channels in a thermal engine Tn that produces a net 
amount of work Wn with probability l — Cn- Let us further 
assume that the probability of failure vanishes in the limit 
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of large n, i.e., lim„_j.oo En = 0. Under these conditions, 
the average asymptotic work W = lim sup satisfies 

n—^oo 

N 

w (8) 

i=l 

As indicated above, this bound is achievable with the 
use of catalysts and a sublinear amount of quantum co¬ 
herence. 

Proof. In any protocol for work extraction, the initial 
state of the system will be given by the catalysts a cat, 
a number of thermal states Tth, and the work system 
in state |0)u,. The initial state of the system is hence 
Po = o-cat®rt;i(8)|0){0|u,, with free energy F{acat)+F{Tth). 

Suppose that now we apply a sequence of energy- 
conserving unitaries. At time t, the state of the overall 
system is pt, and we apply the channel fls(t) over part of 
the whole system, possibly followed by some other ther¬ 
mal operation. Let us analyze how the free energy of pt 
can increase in the above step. Calling Ht the Hamilto¬ 
nian of the whole system, from the definition of IL(r2, H) 
and the additivity of the distillable work we have that: 

(g)l,p0 < ® 1,Ht) = 

(9) 

Now, any intermediate energy-conserving unitary in- 
between the use of any two of the channels will 

keep the free energy of the overall system constant. Call¬ 
ing p the state of the system at the end of the protocol, 
we hence have that 

N 

F{p)<nY,W{n,,H,)+F{acat)+F{Tth). (10) 

From the subadditivity of the von Neumann entropy, it 
follows that F{p) > F{pcat) + F{pth) + F{py,), where 
^cat, Pth, Pw are, respectively, the reduced density matri¬ 
ces of the catalyst, thermal and work systems. 

At the end of the protocol, the catalyst must be re¬ 
generated, i.e., CTcat = (Teat. AlsO, F{pth) > F^Tth)- It 
follows that the free energy of the work system is bounded 

This system is expected to end up in state |1) with 
probability 1 - e„, i.e., Pu, = (1 - e„)|l)(l| -f e„d-. It 
follows that F{pu,) > (1 — e„)IT„ — KBTh{en), with 
h{p) = —pln(p) — (1 — p) ln(l — p). In the asymptotic 
limit, with n —>■ oo, e„ —>■ 0 , the average asymptotic work 
lim sup ^ is hence bounded by 

Note that this bound also holds if the catalysts are 
recovered up to an error, as long as F{acat) — F{acat) < 
o(n). 

□ 


This result allows to quantify the work extraction ca¬ 
pabilities of different channels. One can check, for in¬ 
stance, that no work can be distilled from a dephasing 
channel. Meanwhile, for a two-level system with Hamil¬ 
tonian H = A|I)(I|, E > 0, the channel that takes any 
state to the excited state |I) provides the highest distill¬ 
able work. 

Gaussian channels 

If our target system is infinite dimensional, in princi¬ 
ple there may exist quantum states possessing an infinite 
amount of energy. If we regard such states as unphys¬ 
ical, we should replace the maximization in eq. dSI by 
an optimization over all states of finite energy. The re¬ 
sulting quantity will hence bound the maximum amount 
of work generated in physically conceivable quantum en¬ 
gines, where the overall state of the system always has a 
finite amount of energy. 

In infinite dimensional systems Gaussian quantum 
channels have a special relevance: they are easy to imple¬ 
ment in the lab, and are extensively used to model par¬ 
ticle interactions with a macroscopic environment. They 
are defined as channels which, when composed with the 
identity map, transform Gaussian states into Gaussian 
states, the latter being those states with a Gaussian 
Wigner function [l5| . An m-mode Gaussian state is com¬ 
pletely defined via its displacement vector di = {Ri) and 
covariance matrix 7 ^ = {{Ri — dil,Rj — d_;T}+), where 
= {Qi,Pi,...,QTn,PTn) are the optical 
quadratures. The action of a Gaussian channel is fully 
specified by its action over the displacement vector and 
covariance matrix, given by: 

d^Xd + z,y^XyX^ + Y, (II) 

where Y-\-ia — iX^aX > 0. Here a denotes the symplec- 

tic form a = 1 ^ g ^ . If the Hamiltonian of the 

system under study happens to be a quadratic function of 
the optical quadratures, i.e., H = R^GR-\-h-R, for some 
real symmetric matrix G and real vector h, then the aver¬ 
age energy of a state with displacement vector d and co- 
variance matrix 7 is given by A = i Tr[Gj]-\-(i^G(l-\-h-d. 
States with finite energy hence correspond to states with 
finite first and second moments. 

When the quadratic Hamiltonian has no zero energy 
modes (that is, when G > 0), Proposition [T] allows to 
easily classify generic Gaussian channels according to 
their capacity to generate an infinite amount of work. 
Indeed, for X^GX — G ^ 0, the channel’s distillable 
work is unbounded: this can be seen by inputting a se¬ 
quence of Gaussian states with constant covariance ma¬ 
trix but increasing displacement vector parallel to any 
positive eigenvector of X^GX — G. Gonversely, as we 
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show in the Supplemental Information, for channels sat¬ 
isfying X^GX — G < Q only a finite amount of work can 
be distilled. 

For such channels there is still the dilemma of how 
much work can be extracted. The next Proposition 
greatly simplifies this problem by showing that, for Gaus¬ 
sian channels 17, the maximization in ([3]) can be restricted 
to Gaussian states: 

Proposition 2. Consider a continuous variable quan¬ 
tum system of m modes, with quadratic Hamiltonian H, 
let n be a Gaussian channel mapping m modes to m 
modes, and denote by Q the set of all m-mode Gaussian 
states. Then, 

W{n, H) = max XF{p, Cl). (12) 

P&Q 

The proposition can be proven by combining Lemma [T] 
with the ‘gaussification’ protocol described in [l^, see 
the Supplemental Information. 

Since W{Cl,H) just involves an optimization over a 
finite set of parameters subject to positive semidefinite 
constraints, (in principle) it can be computed exactly for 
any Gaussian channel Cl. 


One application: collapse engines 

In order to address the measurement problem 0 , and, 
independently, the decoherence effects that a quantum 
theory of gravity could impose on the wave-function 18|, 
, different authors have proposed that closed quantum 
systems should evolve according to the Lindblad equation 


with pt the solution of eq. 0 for the initial state po- 

We suppose that the particle under consideration is 
subject to a harmonic potential, i.e., H = F 

and that, despite the GRW dynamics, the bath’s 
temperature T is constant. A physical justification for 
this last assumption is that the temperature-increasing 
GRW dynamics is countered by radiation from the bath 
into outer space. Hence, as a function of time, the tem¬ 
perature will converge to a stationary value T = T^q 
above the temperature of the cosmic microwave back¬ 
ground (GMB) [^. Finally, we suppose that our set of 
resource operations remains the same: that is, in spite 
of the modihed Schrodinger equation 0, we can still 
switch on and off any unitary interaction that commutes 
with the total energy of the system. Notice that the GRW 
dynamics can be modelled by an open system approach, 
where the particle is interacting with some unknown sys¬ 
tem such that the resulting evolution is given by 0. 
From this viewpoint, we are simply assuming that we 
still have the capacity to interact with the system in the 
usual way. 

In these conditions, we wish to find the maximum 
amount of work that a collapse engine could extract if 
it had access to the evolution equation 0 for a finite 
amount of time t. From Proposition this amounts to 
computing lim^^^o ilW{Cl,t,H). 

First, notice that we can (reversibly) evolve the system 
with the Hamiltonians H or —H, since this corresponds 
to a thermal operation. This implies that we can ignore 
the first term in the right hand side of (TT^ . and what 
remains is a Gaussian channel given by 


lp, = -^-[H,p,]-^[X,[X,p,]], (13) 


d ^ d, 7 ^ 7 -f 


0 

0 



(15) 


where X is the position operator for the particle consid¬ 
ered. 

The effect of the non-unitary term is a suppression 
of coherences in the position basis, effectively destroying 
quantum superpositions. The value of the constant A, 
which can be interpreted as the rate at which this local¬ 
ization process occurs, depends on the particular theory 
invoked to justify eg. (1131) . In the Ghirardi-Rhimini- 
Weber (GRW) theory [l7|, the localization process is pos¬ 
tulated to solve the measurement problem in quantum 
mechanics. To achieve this goal and avoid contradictions 
with past experimental results, A must be roughly be¬ 
tween and 10®s“^m“^, according to latest 

estimations (^ . 

Note that the above evolution is non-thermal. Hence, 
it could be used in principle to extract work from noth¬ 
ingness by means of a suitable thermal engine. We will 
call such a hypothetical device a collapse engine. 

To connect this to our previous setting, notice that the 
evolution equation (1131) defines a quantum channel 

Clstipt) = Pt+st, (14) 


It follows that the energy of any input state will in¬ 
crease by AU = \^dt. From Proposition [2J we can 
estimate the entropy increase by just considering Gaus¬ 
sian states. Now, the entropy of a 1-mode Gaussian state 
is an increasing function of the determinant of its covari¬ 
ance matrix [l5| , which, by the above equation, can only 
increase with time. Hence, W{Clst, H) < AU. 

On the other hand, suppose we input a squeezed state 
with 7 = diag(l/r, S^r). Then the determinants of the 
covariance matrices of initial and final states will be hf 
and hf ^ , respectively. The entropy change of the 

state can thus be made as small as desired by increasing 
the value of r, and so the above bound can be saturated, 
leading to W{Clst,H) = AU. Consequently, the maxi¬ 
mum power at which a collapse engine could in principle 
operate is given by 


dW _ h^A 
dt 4m 


(16) 


Using the upper range estimation A ^ 10®s ^m we 
have that a collapse engine powered by a single electron 
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would produce ^ ~ 10“^^ watt. Assuming total control 
over the electrons of a macroscopic sample, one would 
need a kiloton of Hydrogen to power a 40 watt light bulb. 

Conclusion 

In this Letter we addressed the problem of determin¬ 
ing how much work can be extracted from operational 
-as opposed to state- resources. We proved that the so¬ 
lution to this problem in the asymptotic limit is given 
by a single-letter formula that quantifies the amount of 
distillable work that a channel can, in principle, generate 
when supplemented with thermal operations and catalyst 
states. Moreover, we found how this quantity can be de¬ 
termined for bosonic channels, and computed it exactly 
for the case of the GRW dynamics, hence determining 
the maximum power which a hypothetical collapse en¬ 
gine could provide for free. 

Note that we have only studied operational resources 
regarding their capacity to generate work. An interest¬ 
ing topic for future research is to extend our results and 
draw a map of the inter-conversion relations between dif¬ 
ferent operational resources. In the case of state resources 
there is a unique monotonic quantity, the free energy, de¬ 
termining the optimal rates for state transformations Q . 
In this work we have identified an operational monotone, 
the distillable work, but we suspect that there may be 
many others. 
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The deterministic distillable work can be 
superactivated 


At the end of the protocol, the work system will have 
evolved to | 1 ). 

If, rather than implementing step |3l we add the state 
in the preparation stage, then it is trivial to find an 
energy-conserving unitary U 2 that at the last step would 
produce exactly the same amount of work. That is, we 
would have extracted work from the resource state |' 0 ). 

In 0, however, it is shown that no work can be ex¬ 
tracted from a single copy of |' 0 ), even with the use of 
catalysts. This implies that the deterministic distillable 
work of a single use of H is zero. 

Suppose, now, that we have access to two uses of fl. 
Then we can prepare two copies of |'0), from which a non¬ 
zero amount of work can be deterministically extracted 
via thermal operations 

Gaussian channels with finite distillable work 

Let H be a Gaussian channel whose action on the dis¬ 
placement vector d and covariance matrix 7 of the m- 
mode input state is given by 

d^ Xd + z,"/ ^ X-fX^ + Y. (17) 

Suppose that H acts on a system with Hamiltonian H = 
n^GR +h-R, where R is the vector of optical quadra¬ 
tures. Under the assumption that G = G — X"^GX > 0, 
we wonder if the difference between the free energies of 
the input and output states is bounded, i.e., whether 
AF{p, H) < AT, for some K < 00 . 

Call Uq the energy of the input state; and 7 , d, its co- 
variance matrix and displacement vector. Then we have 
that 


Consider the channel H that takes any state of a target 
two-level system with Hamiltonian H = i?i|l)(l| to the 

state |'0) oc ^|0)-f That is, H(p) = IV'XV'I 

for all p. 

Now, any protocol that pretends to extract work out 
via a single use of H can be divided in three steps: 

1. The system is prepared in a state which comprises 
catalysts, thermal states and the work system (in 
state |0)). That is, po = cicat ® nu 0 |0)(0|^,. 


Ao = i Tr[G 7 ] + {d- d^f G{d - do) - E, (18) 

where do = —G~^hl2, and E = S^Gdo- Hence, defining 
Amin > 0 (Amax > 0 ) to be the minimum (maximum) 
eigenvalue of G we have that 


^ (^Tr[7] -P ||d-dof) 

Minax \ ^ / 


Eq -\- E 

A^min 


(19) 


and consequently 


2. We apply an energy-conserving unitary Ui over the 
whole system. 

3. We apply H, hence replacing a subsystem’s state 
by |T). 

4. We apply a second energy-conserving unitary U 2 
over the whole system. 


Tr[ 7 ] <0(Uo),||d]| <0(v^), 

iTr[ 7 ] + ||d-dof > 0 (Uo). ( 20 ) 

We can now bound the energy difference between the 
input and output states. First, note that AE = Ef — Eq 
can be written as: 
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A£; = — Tr 
2 


iG 


{d- do)^G{d- do) + 0{d). ( 21 ) 


Defining Amin > 0 to be the smallest eigenvalue of G, we 
thus arrive at 


< -Amin(^i Tr[7] + ||d - do)f ^ + 0(d) 

<-0{Eo) + 0{^o) = -0{Eo). (22) 

Let us now bound the entropy of the input state: by 
the subadditivity of the von Neumann entropy, S{p) is 
bounded from above by ^(Pi)^ where pi denotes the 
reduced density matrix of each mode i. S'(/ 3 i), in turn, is 
bounded by the von Neumman entropy of the Gaussian 
state with the same first and second moments as pi, i.e., 
a Gaussian state with covariance matrix 7 ^. Note that 
X)™ 1 Tr[ 7 i] = Tr[ 7 ] < 0{Eo), where the last inequality 
is due to eq. (l20|) . In particular, Tr[ 7 i] < 0{Eo) for 
1 = 1 ,..., m. 

The entropy of a 1-mode Gaussian state with covariance 
matrix 7 is given by 


{N + 1) log(7V -hi) - N\og{N) = 0{log{Nj), (23) 

where N = ^ydet{^)/h < Tr[ 7 ]/ 2 fi (this last inequality 

reflects the fact that the geometric mean of 7 ’s two eigen¬ 
values is smaller than their arithmetic mean). It follows 
that 

Sip) < OihgiEo)). (24) 

Putting all together, we have that 

AF(p, n) <AE + KbTS{p) 

<-OiEo)+OilogiEo)). (25) 

The last expression cannot thus take arbitrarily large 
values, and so the distillable work of channel ft is 
bounded. 


Proof of Proposition 2 

Proposition 2 follows straightforwardly from the fol¬ 
lowing Lemma: 

Lemma 2. Let p he an arbitrary state with finite first 
and second moments, and let pa be the unique Gaussian 
state with the same first and second moments. Then, for 
any Gaussian channel D, we have that 


Sipa) - Siftipa)) > Sip) - Siflip)). (26) 


Proof. Let the action of ft over the displacement vec¬ 
tor and covariance matrix be given by eq. dig. Since 
von Neumann entropies remain the same after a unitary 
transformation, without loss of generality we will assume 
that p’s displacement vector is null. Similarly, we will 
take 2 = 0 in the channel description dni of ft. Now, 
let Un be the n-system ‘Gaussification’ transformation 
described in [l^. Galling = Unp^’^Uf,, we have that 


L„ = Y,Sip^)-Siftip,)) 

i^l 

> ^(p^”))-S'(D®”(p(”))) (27) 

= Sip^^) - Sift^^ip^^)) 

= n{Sip) - Siftip))}, (28) 

where the first inequality is due to Lemma 1 in the main 
text, and the equality follows from the fact that, for all 
states CT, 


[/„L!®-(a)C/+ = L!®"(t/„at/t). 


(29) 


This identity follows from three observations: 1) any 
Gaussian channel with 2 : = 0 is the result of applying 
a symplectic unitary Vs over the target system and an 
ancillary Gaussian state uj with zero displacement vec¬ 
tor; 2 ) n copies of w are invariant with respect to a 
Gaussification operation 3) = 0, where 

Ui, = Un ® Un represents the Gaussification of n copies 
of the system target-ancilla. 

Since the von Neumann entropy is continuous in trace 
norm with respect to collections of states with finite sec¬ 
ond moments, by 1 ^, we end up with 


T 1 , 

lim — = lim - S'(pi) - S'(D(pi)) 

n—^oo Ti ' ^ 


n—^00 Tl 




= ,S(pg) - Siftipa)). 


The lemma hence follows from 1)28^ and (15(1)) . 


(30) 

□ 


Now, let p be an arbitrary state with finite energy (and 
thus finite first and second moments), and let pc be the 
Gaussian state with the same first and second moments. 
Then, p and pc have the same average energy, and since 
D is a Gaussian channel the same is true for flip) and 
ftipa). However, by the previous Lemma, the entropic 
change is bigger for pc, and so AFipc,ft) > AFip,ft), 
proving Proposition 2. 
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If we model the bath as a grey body, then it radiates 
energy at a rate of a{T^ — T*), where Tc is the temper¬ 
ature of the CMB and cr is a constant that depends on 
how well isolated the bath is. The power transferred to 
the bath by the GRW dynamics is, on the other hand, 
independent of T and proportional to A. It follows that 
the bath will reach a stationary temperature Te, whose 
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